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Abstract 

We consider the lattice dynamics in the harmonic approximation for a simple hy- 
percubic lattice with arbitrary unit cell. The initial data are random according to a 
probability measure which enforces slow spatial variation on the linear scale e _1 . We 
establish two time regimes. For times of order e -7 , < 7 < 1 , locally the measure 
converges to a Gaussian measure which is space-time stationary with a covariance in- 
herited from the initial (in general, non-Gaussian) measure. For times of order e^ 1 
this local space covariance changes in time and is governed by a semiclassical transport 
equation. 
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weak convergence of measures, semiclassical transport equation. 
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1 Introduction 



For systems consisting of many interacting "particles" as a rule the slow degrees of freedom 
are linked to local conservation laws. For example for a classical fluid, mass, momentum, and 
energy are locally conserved and as a consequence mass, momentum, and energy density are 
the slow degrees of freedom. Thus if the system starts with some general initial conditions, 
one expects the fast degrees of freedom to die out rapidly. Then, in a spatial region which 
on one side contains many particles and on the other side is still small compared to the total 
extent of the system, thus locally, the statistical distribution on phase space is stationary 
under the dynamics within a good approximation. Since the system has not yet reached 
global stationarity, there is still a slow motion of the parameters characterizing the states of 
local stationarity. In our example of a classical fluid, local stationarity coincides with local 
thermal equilibrium and the local equilibrium parameters, density, momentum, and internal 
energy, evolve according to the Euler equations of fluid dynamics. For other systems with 
many particles, in general, it is a difficult task to identify the relevant probability measures 
stationary in time (and usually also in space). 

Such a picture for the dynamics of systems with many particles has theoretical and 
mathematical support. We refer to [TT]. If the dynamics is of Hamiltonian form, the list 
of worked out examples is rather short. One item on the list is lattice dynamics in the 
harmonic approximation, which has been investigated in great detail by R.L. Dobrushin and 
collaborators [2]. We reconsider this model for two reasons. 

(i) The first one is on a conceptual level. In phonon physics it is standard practice to 
use the Wigner function W(t,r,8) as density of phonons with wave number 9 at location 
r and at specified time t . W evolves according to the semiclassical transport equation 



uj(8) being the dispersion relation of the harmonic crystal. As we will establish, W(t, r, 8)5(0— 
6') at fixed r, t encodes the covariance of a Gaussian measure on phase space which is in- 
variant under the lattice dynamics. Thus (jl.lj) can be understood as the equation governing 
the motion of the parameters which characterize the locally stationary measures. We believe 
that in this way the results of Dobrushin et al. become more transparent and, in addition, 
the link to the physics of phonons is provided, see ^2| for a more detailed discussion. 

(ii) The second reason is technically. In the recent years there has been considerable 
progress in understanding the long time limit of the harmonic crystal in infinite volume [3] . 
More precisely one starts with a probability measure fio which is translation invariant and 
has some mixing properties. If fi t denotes the time-evolved measure at time t , then the 
limit 



is established, where is a suitable Gaussian measure with mean zero. It turns out 
that the techniques for proving (|1.2|) transcribe to the locally stationary situation. Thereby 
the conditions in the work of Dobrushin et al. are considerably streamlined and the proof 
is simplified. We also generalize from one to an arbitrary space dimension and from one 
particle per unit cell to an arbitrary number. 
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In a recent paper jHj, A. Mielke studies the same model and also obtains the semiclassical 
transport equation for the Wigner function. However, Mielke imposes deterministic 

initial data of slow variation, while we impose random initial data with rather strong mixing 
properties. Therefore the results are disjoint and so are the techniques for proving them. It 
is of interest to understand whether a "supertheorem" encompassing both cases has a chance 
to be valid. 

2 Lattice dynamics in the harmonic approximation 
2.1 The model 

We consider a Bravais lattice with a unit cell which contains a finite number of atoms. For 
notational simplicity the Bravais lattice is assumed to be simple hypercubic. Let x € Z d 
and let u(x) be the field of displacements in cell x from the equilibrium position. If u is 
small, we may expand the forces to linear order, which then yields the linearn— component 
discrete wave equation 

u(x, t) = - V(x - y)u(y, t), u(x)\ t=0 = u {x), u{x)\ t=0 = v (x), x E 7L d . (2.1) 

Here u(x,t) = (ui(x, t), . . . , u n (x, £)), u = (u i(x), . . . , u 0n (x)) G M™ and correspondingly 
for Vq(x) . Physically n = dx (number of atoms in the unit cell). Here we take n to be an 
arbitrary positive integer. V(x) is an n x n matrix. The dynamics (|2.1j) is invariant under 
lattice translations. 

Let us denote by Y(t) = (Y°(t),Y\t)) = («(•,*),«(•,*)), Y = (Y °,Y l ) = W')^(-)) ■ 
Then (j2.1j) takes the form of an evolution equation 

Y(t)=AY(t), teM, Y(0)=Y . (2.2) 

Formally, this is a linear Hamiltonian system, since 

AY = j(j Q l^jY=JVH(Y), ^=(_° 1 J), (2-3) 

with the Hamiltonian functional 

H(Y) = ^(v,v) + ^(Vu,u), Y = (u,v), (2.4) 

where V is the convolution operator with the matrix kernel V , the kinetic energy is given by 
—{v, v) — — \ v ( x )\ 2 ) an d the potential energy by -(Vu,u) = — u{x) ■ V(x — y)u(y) . 

x£Z d ~ x,yeZ d 

Here " • " stands for the scalar product in the Euclidean space M. n , resp. in M d . 

We assume that the initial datum Y belongs to the phase space 7i a for some a£E. 

Definition 2.1 7i a is the Hilbert space of pairs Y = (u,v) of M. n -valued functions on Z d 
equipped with the norm 

\\Y\\l = (K^| 2 + K^)| 2 )(1 + M 2 ) Q <oo. (2.5) 

x&Z d 

TC a is equipped with the Borel a -algebra B(H a ) ■ 
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We impose the following conditions on the matrix V. 

El There exist constants C, a > such that \\V{z)\\ < Ce~ a ^ for z G Z d , \\V{z)\\ 
denoting the matrix norm. 

Let V(9) be the Fourier transform of V(x) , with the convention 

V(9) = ^V(z)e iz - e , 6 (ET d , (2.6) 

z&Z d 

T d the d -torus R d /(2irZ) d . 

E2 V is even, in the sense that V(—z) = V(z)* G K, for z G r L d , where V* denotes the 
adjoint of the matrix V as acting on C n . 

Both conditions imply that V(9) is a real- analytic Hermitian matrix- valued function in 
9 G T d . 

E3 The matrix V{9) is non-negative definite for every 9 G T d . 
Let us define the Hermitian non-negative definite matrix 

Q(0) = {V(6)) 1/2 > 0. (2.7) 

fi(#) has the eigenvalues < uji{6) < u 2 (9) . . . < oj s {6) , s < n and the corresponding 
spectral projections H a (0) with multiplicity r a = tril o -(0) . 9 i— > co , cr (6') is the cr-th band 
function. There are special points in T d , where the bands cross, which means that s and 
rv jump to some other value. Away from such crossing points s and r a are independent of 
9 . More precisely one has the following lemma. 

Lemma 2.2 (see J2J Lemma 2.2]). Let the conditions El, E2 hold. Then there exists a 
closed subset C* C T d such that 

i) the Lebesgue measure of C is zero. 

ii) For every point G T d \ there exists a neighborhood 0(Q) such that each band 
function ui a {9) can be chosen as real- analytic function in 0{Q) . 

Hi) The eigenvalue iO a {9) has constant multiplicity in T d \ . 
iv) For 9 G T d \ C* , the spectral decomposition 

s 

n(9) = 5> CT (0)n CT (0) (2.8) 

holds, where n CT (^) is an orthogonal projection in M. n . n CT is a real- analytic function on 
T d \C*. 

For 9 G T d \ C* we denote by Hess {oJ a ) the matrix of second partial derivatives. Our 
next condition is the following. 

E4 Let D a {9) = det (B.ess(u a (9))) . Then D a does not vanish identically on T d \ , 

(T=1,...,S. 

Let 

C = {9 G T d : det V(9) = 0} and C a = {9 eT d \C* : det ( Hess = 0}, a = 1, . . . , s. 

(2.9) 

The following lemma has been proved in [3J Appendix]. 
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Lemma 2.3 Let the conditions El - E4 hold. Then the Lebesgue measure of C k vanishes, 
k = 0, 1, s. 

Our final conditions on V are the following: 

E5 For each a ^ a' , u a ± ou a > does not take a constant value on 9 G T d \ . 
This condition holds trivially in case n = 1 . 
E6 H^- 1 ^)!! G L l {T d ). 

If C = , then || t> 1 (6>) || is bounded and E6 holds trivially. 

Remark 2.4 The conditions El - E6 are fairly general. In particular they can be checked 
for the case of nearest neighbor coupling only, for which 

n d 

(Vu,u) = ^2 ^2 (5^7*|wfc(^ + e<) - u k (x)\ 2 + m 2 k \u k (x)\ 2 ^ , ^ k > 0, m k > , (2.10) 

k=l x& d i=l 

where = (8n, . . . , 5^) • Then the eigenvalues of V{9) are 

u k (9) = -^27! (1 - cos 0i ) + ... + 2 7d (l - cos9 d ) + m\ . (2.11) 

These eigenvalues still have to be labelled according to magnitude and degeneracy as in 
Lemma l2~2l Clearly El - E5 hold. In case all m k > the set Co is empty and condition E6 
holds automatically. Otherwise, if m k = for some k , Co = {0} . Then E6 is equivalent 
to the condition (9) G L 1 (T d ) , which holds if d > 3 . Therefore, the conditions El - E6 
hold for ()2.10j) provided either i) d > 3 , or ii) d— 1,2 and all m k > . 

Proposition 2.5 (see J3J Proposition 2.5]). Let El and E2 hold and choose some a G R. 
T/ien 

/or am/ Y"o G i/jere exists a unique solution Y(t) G C(R, 7i a ) to £/ie Cauchy problem 
T/ie operator U{t) : Y i— > zs continuous in 7i a , \\U{t)Yo\\ a < C{t)\\Yo\\ a . 

2.2 Random initial data 

We assume that Yo is a random function with distribution fio ■ 

Definition 2.6 fit is a Borel probability measure in 7i a which gives the distribution of 
Y(t) , 

IH(B) = fi (U(-t)B), BeB(H a ),teR. 
Expectation with respect to fit is denoted by E t . 

We set V = D@D with D = C (Z d ) ®M n , where C {Z d ) denotes a space of real sequences 
with finite support, and (Y,V) = (Y°,^°) + (Y 1 ,^ 1 ) for Y = (Y°, Y l ) G H a and * = 
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(\I/ , SI/ 1 ) G T> . For a probability measure \i on 7i a we denote by fi the characteristic 
functional (Fourier transform), 



= y exp(z(Y,*)) /i(c/y), *GP. 



A measure /i is called Gaussian of zero mean, if its characteristic functional has the form 

/2(*)=exp[-ifi(*,*)], *eP, 

where Q is a real non-negative quadratic form on D . A measure /i is called translation 
invariant if fi(T h B) = fi(B), B G B(H a ), h G Z d , where T h F(x) = Y(x - h) , x G Z d . 

Let O(r) denote the set of all pairs of subsets A, B C Z d at a distance dist (A, B) > 
r and let a (A) be the a -algebra in H a generated by Y(x) with x £ A. Define the 
Ibragimov-Linnik mixing coefficient of a probability measure fi on 7i a by (cf. [HI Definition 
17.2.2]) 

lu(AnS) -u(A)u(B)\ 
<p(r)= sup sup ^ J —p ^ K l X (2.12) 

(A,B)eO(r) Aaa{A),B aa{B) 
H(B) > 

Definition 2.7 ^4 measure \x satisfies the strong uniform Ibragimov-Linnik mixing condi- 
tion if tp(r) — > as r -> oo . 

3 Main results 

3.1 Spatially homogeneous initial measure 

In this subsection we assume that the initial measure, hq , is spatially translation invariant 
with the following properties. 

51 Yq(x) has zero expectation value, 

E (y (x))=o, xez d . 

52 fi has translation invariant correlation matrices, i.e., for i, j = 0, 1 , 

Qi j (x,x') = E (Y*(x) ® Yj(x')) = q${x - x'), x,x' G Z d . (3.1) 

Here for a, b, c G C n we denote by a <8> 6 the linear operator (a <g) 6)c = a Y^=i ^j c j ■ 

53 fiQ has a finite variance and finite mean energy density, 

e = E (|F °(x)| 2 + IFq 1 ^)! 2 ) = tr g 00 (0) + trg o n (0) < oo, x G Z d . 

54 fi satisfies the strong uniform Ibragimov-Linnik mixing condition with 



CO 



d - l ip l,2 {r)dr< oo. (3.2) 



o 



In [3] we prove the weak convergence of the measures fit to a limit measure /ioo on the 
Hilbert space TL a with a < —d/2, which means 

lim J f(Y)m(dY) = J fWn^dY) (3.3) 

for all bounded continuous functions / on TC a . //qo is a Gaussian measure on 7i a . 

Theorem 3.1 (see J3J/J. Let d,n> 1 , a < —d/2 , and assume that the conditions El - E6 
and SI - S4 hold. Then 

i) the correlation matrices of the measures /i f converge to a limit, for i,j = 0, 1 , 

Q^(x,x') = J (Y\x) ®Yi(x')) fi t (dY) ^ Q^(x,x'), t - oo. 

ii) The convergence in hS. holds. 

Hi) The limit measure fj,^ is a Gaussian measure on 7i a ■ 

iv) The correlation matrix of /Jm is translation invariant, Qoo(x,x') = qoo(x — x') , and has 
the Fourier transform 



(7=1 



where H a (9) is the spectral projection from Lemma \2.^\ iv) and 



M o (0) = \ (»(») + C(S)q«(S)C(S)') 

0(0) - ( ° "'I" ) • (3.4) 



with 

n(a\ — I 

-Q(9) 

v) The measure /ioo is time stationary, i.e. [U(t)]*[/,oo = fioo , t£R. 

The projection of the initial covariance %{9) to the limiting covariance can be 

stated more concisely through introducing the complex-valued field 

a(x) = (v 1,A u{x) + iV~ 1/4 v(x)^ G C n , xeZ d , 

with complex conjugate field a(x)* and distributional Fourier transform d(9) . 

Obviously E t (a(x)) =0. The covariance has two parts. By Theorem 13.11 the aa-, 
equivalently the a*a*-, covariance satisfies 

lim E t (a(x) ® a(x')) = 0. 

For the a* a -covariance we define 

Et (a(ey ® a{9')) = (2ir) d 5{9 - 9')W(t, 9) , 
using the translation invariance of [i t . Note that W(t, 9) > . Then 

s 

lim W{t, 6) = Y* UJ9)W(0, 9)U a (9) . (3.5) 



(7=1 
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3.2 Initial measure with slow variation 

Let {/j,q,£ > 0} be a family of initial measures. Roughly, in a linear region of size e~ x , 
£<1, /Xq looks like the spatially homogeneous initial measure from Section 13.11 However 
the covariance Qq depends on the spatial region under consideration, and not only on the 
difference x — x' . 

To be more precise let us introduce the complex In x In matrix-valued function R on 
R d x T d , through 

with the following properties. 

II For every fixed r6l d and i, j — 0, 1, the entries of the matrix-valued function R are 
bounded on T d and the inverse Fourier transform 



R l] (r,x) = (2tt) 
satisfies the bound 

\R ij {r,x)\<C{l + \x\)-\ xeZ d , (3.6) 
where C is some positive constant, 7 > d. 

12 For every fixed r GK d , the matrix-valued function R satisfies 

R 00 {r,e)>0, R 11 {r,6)>0 (3.7) 

R ol {r,0) = R lo (r,6)*, 9eT d . 

13 For every fixed r e R d and e T d , the matrix i?(r, #) is non-negative definite. 

14 For every 8 ET d , R^(-, 6) , z, j = 0, 1 , are C 1 functions and the function 

sup max (\R ij (r,6)\, I V r R ij (r, 0)1) 

j £ Td»,i=0,l 

is bounded uniformly on bounded sets. 

Let Eq stand for expectation w.r.t. the measure /xg • We assume that 

E £ o (^(a;))=0 (3.8) 

and define the covariance 

Qi j (x,x') =E £ (Y i (x) ®Y j (x')), x,x'eZ d , i,j = 0,1. 



r — > 

eel 
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Definition 3.2 We call a family of measures {/Zq, e > 0} a family of slow variation for R 
if {Q l J(x,x'),e > 0} satisfies the conditions VI - V4 listed below. 

VI For any e > there exists an even integer N £ such that 

i) for all I6l d and x, x' G I M , 

\Qi j (x,x')-R ij (eM,x-x')\ < Cmin[(l + \x - x'|)~ 7 , eN e ], (3.9) 

where C , 7 are the constants from (|3.6|) . and /a/ is the cube centered at the point M with 
edge length N e , 

I M = {x = ( Xl , ...,x d )eZ d : \ Xj - Mj\ < N £ /2, M = {M u . . .,M d )}. (3.10) 

ii) N £ ~ e~P as e -> , with some /3 G (1/2, 1) . 
V2 For any e > and all x, x' G Z d , z, j = 0, 1 , 

\Qi j (x,x')\ <C(l + \x-x'\)-i 

with constants C , 7 as in (|3.6j) . 

V3 For any e > and any ^i,^ G V with dist (supp supp > P > there exist 
constants C > and ft G (0, 1) such that 

|Eg( e i<y,*i) e i(y,* a >) -Eg( e W0)Eg( e Wa>)| <c(i + P )- K . 
V4 For correlation functions of the fourth order 

Mi 4) (x\x 2 ,x 3 ,x 4 ) = Eg(y(x x ) ® F(x 2 ) ® F(x 3 ) (g>F(x 4 )), a; 1 , . . . , x 4 G Z d , 
we require that 

|M £ (4) (x\aV^ 4 )l < C Yl (l + \x h -x i2 \)^(l + \x i3 -x u \y\ 

(i 1 ,i 2 ,t3,*4)eP{l,2,3,4} 

where P{1, 2, 3, 4} is a permutation of the numbers 1, 2, 3, 4 , and 7 > d . 

Definition 3.3 /if is a Borel probability measure in 7i a which gives the joint distribution 
ofY(t), 

(4(B) = fi £ (U(-t)B), B G B(H a ), teR. 
m/ 1 T/ie correlation functions of the measure \i\ are defined by 

Q%(x,y) = jY\x)®Yi(y) f il(dY)=W (Y i (x,t)®YJ( y ,t)), ij = 0,1, x,yeZ d . 

Here Y l (x,t) are the components of the solution Y(t) = (Y°(-, t), Y l {-, t)) . 
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3.3 Covariance in the kinetic scaling limit 

The family fi £ Q , s > , of initial measures has slow spatial variation on scale e -1 and for 
long times, roughly of order e~" f , < 7 < 1 , in essence Theorem 13.11 applies locally, which 
implies that locally the projected measure is attained. This measure is then almost invariant 
under the time evolution. Thus one needs a time span of order rje , r 7^ , to see changes 
in the projected part of the covariance. 

To state a precise result we introduce the scaled n x n Wigner matrix through 

W%t; r,9) = Yl ^ E r/e( a *([ £ ' lr + Vl 2 \) ® - I//2])) ■ (3.11) 

y&Z d 

By our assumptions on //q , the following limit exists 

\imW £ (0;r,9) = - (V/ 2 £ 00 (r, 9)Q}' 2 + tl^R 11 ^, 9)n~ 1 ' 2 

+iQ 1 / 2 R°\r, 9)Q-^ 2 - iQ-^ 2 R 10 (r, 9)n 1 ' 2 ^ 
= W(O;r,0). 

We also define the projected initial Wigner matrix, compare with (|3.5)L 

s 

W p (r, 9) = J2 n - (0)W(0; r, 9)U a (9) (3.12) 
0=1 

and its time evolution 

s 

WV(t; r,9) = J2 ^(0)W(0; r - tVu^O), 9)U a (9) . (3.13) 

<T=1 

Theorem 3.4 Let the conditions VI - V2 and El - E6 hold. Then for any r e M. d and 

r/0 the following limit exists in the sense of distributions, 

\imW e (T;r } 9) = W p (T;r,9) . (3.14) 
In addition, for the remaining part of the covariance, 

lim^V^E^aQ^V + y/2])®a([£- 1 r -y/2])) =0. (3.15) 

y&Z d 

We remark that in the a -th band the Wigner function evolves according to the transport 
equation 

J^(r, 9) + Vu a {9) ■ V r / t (r, 9) = , (3.16) 

where the initial conditions are given by the initial Wigner matrix projected onto the a -th 
band. 

The conditions VI and V2 on the initial measure are written in position space. Therefore 
it is natural to prove the limiting covariance first in position space, which will be stated in 
Theorem 14.11 From it we deduce the limiting Wigner function of Theorem 13.41 
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3.4 Local stationarity 

So far we studied only the covariance. A more detailed statistical information is provided 
by considering the random field Y at the kinetic time r/e , r^O, and close to the spatial 
point [r/e] G Z d . For this purpose let T h , h e Z d , be the group of space translations. The 
measure at r/e is then defined through 

f4/e,r = T -[r/e]^T/e • ( 3 -17) 



Theorem 3.5 Let the conditions VI - V4 and El - E6 hold. Then for r^O, in the sense 
of weak convergence on Ti. Q , 

!im/4/ e,r = ^T,V ( 3 - 18 ) 

fi^r is a Gaussian measure on Tt a , which is invariant under the space translations Th and 
time translation U (t) . jj^ r has mean zero and covariance 

q%{x - x') = E^Y^x) ® Y j (x')) , 
expectation with respect to y^ r . The covariance is determined through W p (r; r, 6) as 

mCW) = nW'^W = \(w*{r-r,9) + w*(r, -ey) (3.19) 

and 

t]M = -C( ) = -^(W p (r;r,9) - W»(T;r,-e)*) . (3.20) 

We conclude that close to r/e in space and close to r/e in time the random field Y^(x, t) 
is a stationary Gaussian field. Its distribution at fixed local time t is given by fi^ r while in 
time it evolves deterministically according U (t) . In this sense locally in space and time the 
random field is stationary with statistics determined through the Wigner matrix at (r, r) 
and the microscopic dynamics, compare with ()3.19|h ([3.20)1 . 

4 Convergence of correlation functions 

At first we introduce the matrix q T . r {x) . In Fourier space, 

s 

<T=\ 

where ITo-^) is the spectral projection introduced in Lemma f2. 21 iv) , 

M£(r;r,0) = §(R^(<r;r,0) + C(9)R° + (T;r,9)C*(9)), 
M*(r;r,0) = ±(C(0)R^(t; r, 9) - R*(r; r, 9)C*(9)), 

with matrix C{9) as in (|3.4j) and 

R±(r;r,0) = ~(R(r + Vu a (9)r,9) ± R(r - Vw ff (0)r,0)). 
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(4.1) 

(4.2) 
(4.3) 



Theorem 4.1 Let the conditions VI - V2 and El - E6 hold. Then for any r E M. d , 

x, y G Z d , r^O £/ie correlation functions of measures r converge to a limit, 

SS^r/ e (t r /e] + ac, [r/e] + y) = g^(s - y). (4.4) 

We outline the strategy of the proof. For the proof we use the cutting strategy from 
|3| combined with some techniques from |2] , where Theorem 14.11 has proved for the case 
d = n — 1 (see [21 Theorem 3.1]). Note that in j2] it is assumed the stronger conditions on 
matrix V than E3, E4, namely, u>(9) > , and the set 

{6e [-tt.tt] : u"(9) = u"'{9) = 0} 

is empty. Under these conditions, in [2] the uniform asymptotics of the Green function is 
proved, 

sup|£ t (a;)|<C(l + |i|)- 1/3 . (4.5) 

x£Z d 

This bound plays an important role in the proof of However, if n > 1 , then uo s 
may be non-smooth because of band crossing, and if d > 1 , the set where the Hessian 
vanishes does not consist of isolated points. Therefore a strong estimate as ()4.5|) is unlikely 
to be valid, in general. To cope with such a situation, we split Gt(x) into two summands: 
Qt{x) = G(( x ) +Gt( x ) > where G({x) has a support in the neighborhood of a a "critical set" 
C C T d , and Gf( x ) vanishes in the neighborhood of C . The set C includes all points 9 G T d 
either with a degenerate Hessian of u a (9) , or with non-smooth uj a {9) (see Definition 14.11)) . 
We show that the contribution of Q{ \x) is negligible uniformly in t (see (14.15(0 . Hence, it 
allows us to represent correlations functions Q £ ,t/s in the form: Q £jT / £ = Q g £ T / £ + Q r £ T j e , 
such that 

x',y'<=:l, d 

For the remainder Q r £T / £ = Q £ , T /e ~ Q 9 eT /e we P rove that Q r £T / £ (x,y) = o(l) uniformly in 
r ^ 0, e > and x, y G Z d . The last fact follows from two key observations: i) mesC = 
(Lemma I2.2J1 and ii) the correlation quadratic form is continuous in £ 2 , see Corollary 14.31 
Up to this point we apply the "cutting strategy" from jSlHj- Finally, in Section 4.3 we prove 
that Qg T / e ([r/s] +x,[r/e)+y) converges to a limit as e — > , using the techniques of 0. In 

addition, the asymptotics of Gf {%) > ( see Lemma I43|) of the form Gf{ x ) ~ (1 + \ t\)^ d ^ 2 plays 
the important role, since it replaces the asymptotics ()4.5|) and also simplifies some steps of 
the proof of [2|. 

4.1 Bounds for initial covariance 

Definition 4.2 By £ p = £ p (Z d ) ®M n , p > 1 , n > 1 , we denote the space of sequences 

( \ 1//p 

f{x) = (fi(x),...,f n (x)) endowed with norm \\f\\ p = f \f( x )\ P ) 
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Lemma 4.3 Let condition V2 hold. Then for i,j = 0, 1, the following bounds hold 

\Qe( x >v)\ ^ C <oo for all iG Z d , 

\Qe( x >y)\ ^ C <oo for all ye Z d . 



Here the constant C does not depend on x, y EZ d and e > . 

Corollary 4.4 Lemma \^~S\ implies, by the Shur lemma, that for any $, $ G £ 2 the following 
bound holds: 

\(Q e (x,y),${x)®*(y))\ <C||$|| 2 ||tt|| 2 , 
where a constant C does not depend on e > . 

4.2 Stationary phase method 

Applying Fourier transform to ()2.2j) we obtain 

Y(t) = A(G)Y(t), t G R, Y(0) = Yq. (4.6) 

Here we denote 

The solution to (|2.2|) admits the representation 



Y(x,t) = Y,St{x-y)Y (y), (U 



where the Green function Qt(x) has the form 

Gt(x) = {2nY d J e^exp (A(6)t) d6. 

Note that 

n i f\\ ( cos fit sin fit Vr l \ , . 

g ^={-smntn cos fit J' ( 4 - 9 ) 

where fi = Q(6) is the Hermitian matrix defined by (|2.7p . Hence, we can rewrite Qt(x) in 
the form 

g t (x)=^2 f e- ix - d e ±iuJ ° {d)t a±(8)d6. (4.10) 
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We are going to apply the stationary phase arguments to the integral (j4.10j) which require a 
smoothness in 9 . Then we have to choose certain smooth branches of the functions a^{9) 
and oo a {9) and cut off all singularities. First, we introduce the critical set as 

c = c,LKUU{«^%?M. («i> 

(7 = 1 7 = 1 (7 = 1 1 

with C* as in Lemma 12721 and sets Co and C a defined by ()2.9j) . Obviously mesC = 0. 
Secondly, fix an 5 > and choose a finite partition of unity 

M 

f(9)+g(9) = l, g(9) = Y,9m(0), OeT 1 , 

771=1 

where /, g m are non- negative functions from (T d ) , and 

supp / C {9 G T d : dist(0, C) < 5}, supp g m C {9 <E T d : dist(0, C) > 5/2}. (4.12) 
Then we represent Qt(x) in the form Qt(x) = Q({x) + Gti x ) > where 

<?/(*) = (2n)- d J e~ lx - e f(9)g t (9)d9, (4.13) 

e 4 **g(0)§ t {O)d0 = EE/ ^(0)e-^ e± ^ (e) ^(0)^.(4.14) 

By Lemma 12.21 and the compactness arguments, we can choose the supports of g m so small 
that the eigenvalues u a (9) and the amplitudes ag(9) are real-analytic functions inside the 
supp g m for every m . (We do not label the functions by the index m to not overburden the 
notations.) For the function G{ (x) , the Parseval identity, (|4.9jl . and condition E6 imply 

\\$l{-)\\l = C J \Gt{9)\ 2 \f{Q)?d9<C J |£ t (0)| 2 d0^O as 5^0, (4.15) 

T d dist(0,C)<5 

uniformly in tel. For the function Gf(x) the following lemma holds. 
Lemma 4.5 Let conditions El - E4 and E6 hold. Then 
i) sup \G 9 t {x)\ < C r d/2 . (4.16) 

x&L d 

ii) For any p > t/iere e:nsi C p ,-y g > swc/j £/iat 

|Sf(rr)| < C p (\t\ + \x\ + l)- p , |x| > 7fl t. (4.17) 

Proof Consider Qf(x) along each ray x = vt with arbitrary v G M. d . By f!4.14|) . one obtains 

M s 



$ 9 M) = E E / 

771=1 ±,(T=l r „j 



-*(*-*tM*))* o ±(0) 
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This is a sum of oscillatory integrals with the phase functions 4>^(0) = 9 ■ v =F ^<r(#) • Since 
uJcr(8) is real- analytic, each function 0^ has no more than a finite number of stationary 
points 9 G suppg, which are solutions to the equation v = ±Vuj a (9) . The stationary 
points are non-degenerate for 9 G supp g m , by (j4.11|) and (|4.12|) . since 

det (^w) = TDa{6) ^ °' 9 G supp ^- 

Therefore, Q 9 t {vt) = 0(t~ d l 2 ) according to the standard stationary phase method [7| ITU]. 
This implies the bounds (j4.16|) in each cone \x\ < ct with any finite c. 

Further, denote by v g = max max max |Vuv(#)| . Then for |u| > v g the stationary 

m <x=l,...,s SSsupp g m 

points do not exist on the supp g . Hence, the integration by parts as in flO] yields Qf(vt) = 
0(t~ p ) for any p > 0. On the other hand, the integration by parts (see ()4.13jl ) implies 
the similar bound Gf(x) = 0{[t/\x\) l ^j for any I > 0. Therefore, (j4.17j) follows with any 
7 9 > v g , which means that the bounds (|4.16|) hold everywhere. ■ 

4.3 Proof of Theorem 14.11 

Proof Step i). The representation (j4.8j) gives 

Q £>t (x,y)=El(Y(x,t)®Y(y,t)) = £ g t (x-x')Q e (x' ,y')g t (y-y')* (4.18) 

x',y'GZ d 

for any t G R 1 . Corollary 14.41 and ()4.15j) imply 

Qe,t( x >y)= E < 5t^ x ~ x ')Qei. x \y')Gt{y-y'T + o(i), 

x',y'&1 d 

where o(l) — > as 5 — > uniformly in t G 1 and x,y E Z d . In particular, setting t = r/e , 
x = [r/e] + I and y = [r/e] we get 

g., T /«([r/ e ]+Z,[r/ e ]) = £ Gl /£ {[r/e]+l-x')Q £ {x\y')gi /£ {{r/e]-yr + o{l) 

x',y'eZ d 

= E + x ')Qe([r/s]-x', [r/e]-y')G 9 r/£ (y'y + o(l). 

x',y'GZ d 

Let c = 7 g + \l\ . Then Lemma f4. 51 ii) and condition V2 imply 

g. |T /e([r/e]+Z, [r/e]) = £ # /e (J +*')Q £ ([r/e]-x', [r/e]-^)^^)* 

x',y'e[-cr/e,CT/e] d nZ d 

+n(e,r) + o(l), 

where lime _:P ri(£:, r) = for any p > and r G M 1 . 

S*tep mJ. We divide the cube [—CT/e,cr/e] d onto the cubes I n N E (see (I3.1U|) ). 

[-CT/E,cr/E] d C [J 4jv £ , 

neJ 
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where J = {n = (m, . . . , n</) G %> d , \nj\ < [cr/(eN e )] + 1} . Then 

Q e , T/E ([r/s]+l,[r/s\) = £ £ 0* /e (Z +x')Q £ ([r/e] [r/e] -j/)^ /e {r/) m 

m,neJ x ' e / mJVe 
2/' G I„jv e 

+ri(e,r) + o(l) 
me J x',y'el m N e 

+r 2 (e,r) + n(e,r) + o(l), 

where 

r 2 (e,r) = £ ^ /£ (/+a;')g £ ([r/e]-a;',[r/,]-y')^ /£ (y / )*- (4-19) 

m, n S J, m ^ n 
x' e / m ]v e , y' e l„N e 

Now we prove that 

r 2 (e,r)^0 as e^O (4.20) 

for any r G M 1 . Indeed, we divide the sum in the RHS of ()4.19|) onto two sums S\ and 
S2 , where the first sum Si is taken over all x' G I m N £ and y' G J n jv £ and m,n G J such 
that 3j G {1, . . . , d} : \rrij — rij\ > 2 ; the sum S2 is taken over all x' G I m N e and y' G I„jv e 
and m, ra G J such that m ^ n and Vj = 1, . . . , d : |mj — rij\ < 1 . By Lemma f4.5[ i) and 
condition V2, the sum S"i is estimated by 

C(l + T/e)- d (r/e) d £ (l + N)" 7 , 

s&Z d ,\s\>N £ 

which vanishes as e — > , since N £ — > +00 and 7 > d . To estimate the second sum (the 
contribution of nearest neighbors I m N e and I n ^ e ) we choose a number p > d+ 1 and divide 
the sum onto two sums: S 2 = £21 + 5*22 , where the sum S21 is taken over all m G J and 
x' G ImN e 1 n G G J : n 7^ m, Vj : |mj — nj| < 1} and y' G J„Ar e such that \x' — y'\ > N^ p 
and the second sum S 2 2 is taken, respectively, over y' such that W -y'\ < Ne /p . The 
contribution of "non-boundary zones" S21 is 

C(l + r/e)- d (r/e) d £ (1 + \s\)^ 

s& d ,\s\>Nc /p 

which vanishes as e — > . The contribution of "boundary zones" S22 is order of 

C(l + T/e)- d {T/eN £ ) d N l J p+d - 1 N d/p ~ CW^ 1 ^ -1 . (4.21) 

The number p is chosen such that (d + l)/p — 1 < . Hence, f)4.21|) vanishes as e — > by 
condition VI, ii). The decay ()4.2()j) is proved. 

Step in). Now we can apply the condition VI, i) at the points [r/e] — x', [r/e] — y' of 
the same cube I\ r / E ]-mN e and obtain 

\Q e {[r/e]-x',[r/e]-y') -R{e[r/e] - emN e ,y' - x')\ < Cmin[(l + [x' - y'\Y\ eN £ ]. 
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Then 

Q £>T/e ([r/e)+l,[r/e)) = £ ^ Q 9 r/e {l +^')R{e[r /e] - emN^y' - x')G a T/e {y'y 

+r 3 (e, r) + r 2 (e, r) + n(e, r) + o(l). (4.22) 
Let us prove that lim r 3 (£,r) = for any r G M 1 . Indeed, since for fixed x' G I m N e the 

e— >0 

sum min[(l + |x' — y'|)~ 7 , £-/V e ] is the order of {eN £ ) l ~ d ^ , we get, by Lemma f4. 51 i), 

\rs(e,r)\ < £ 1^ +*')| niin[(l + \x' - y'\Y\ eN £ ) \g 9 T/e {y')*\ 

m&J x>,y'£l mNe 

< C(l + T/e)- d {T/(eN e )) d N d £ (eN e ) 1 - dh ~ e (i-fl(i-«i/<y) _ , e _, 0, 

by condition VI ii), since /3 < 1 and 7 > d . 

Step iv). By similar arguments, as in steps i) and ii), the sums in the RHS of (|4.22j) can 
be taken over {y 1 G Z d , m G J, x' G / m Ar e } • The sum in y' is a convolution which can be 
expressed by the product in the Fourier transform: 

Q £ ,r/e([r/e)+l,[r/e}) = (2^ £ £ je-^e^< d '-^g 9 T/e {d)R{e[r/e]-emN £ ,d') 

xg 9 T/£ (6')*d6d6' + o T (l) + o(l), 

where o T (l) — > as £ — > for any r G IR 1 \ {0}. Further, since 7 m 7v e = {x' G Z d : 
(rrij - l/2)N e < x) < {rrij + l/2)N £ , j = l,...,d}, then 



x'eImN s j = 1 



where F{9 j ,N £ ,m j ) = e ^N e ( mj +i/2) _ ^6^^-1/2) _ changing variables i) (0,0') -> (2, 0') , 
z = 0' -6 and ii) (2, 0') (s, 0) , = 0', one obtains 

Q e , T /e([r/ e ]+Ur/e]) = M" 2 ^ / e^-'fl ^fegglj^ ^y g _ ,) 

xB(e[r/e] - emiV £ , e)Q a T/e {9)* d9dz + o T (l), (4.23) 

where a(z) = — if z G (— 7r, 7r) \ and a(0) = 1 . Note that 

Ql(Q) = g(0)(costt{6)t + sinQ{6)t C{6)) . 
Hence, in the integrand in ()4.23|) we have for t = t/e , 
G 9 t {0 - z)R(e[r/e] - emN £ , 0)0?(0)* 

s 

= n*(0-z)g(e-z)(cosu a (6-z)t + smu a (6-z)tC a (6-z)) 

R(e[r/e] - emN E , 6)g{6) (cos uv(0)t + sinuv(0)£ C*,(0))lLy(0), (4.24) 
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where C a (9) = ( ^ /n\ ^ ) . Let us consider one of the terms in (|4.23j) . The 



-u>M o 

proof for the remaining terms is similar 



7T 

i 1 f ^-L^(9)r/ £ „/m V" ( [ e iz d i d a ( z d)F(z d ,N £1 m d ) 



I £ = (2vr)- M i I e- m e^'^g(9)J2(f 



2 / 2 a(z 2 )F(z 2 ,N £ ,m 2 ) / /~ ^ otfcQFfci, iV E , mQ c±Mr( g_ z)r / £ 

^2 ^ i ^1 



-7T 



xp(0-z)n ff (0-^)^(e[r/e] - £mAT £ , 0)n CT /(0)dziJrfz 2 J ...dz d )d9. (4.25) 

Introduce 2/1 = v 1 {9 u B 2 - z 2 , . ..) = ±[ViW <T (0i, 2 - z 2 , ■ ■ • )r/(eJV s )] , ^2 = ^2(^1,^2,^3 - 
z 3 , ■ • • ) = ±^2^(0!, 3j 3 - 23, -WO^)] , -, ^ = ^(0) = ±[V d w <T (0)r/(eiV e )] . 



Lemma 4.6 Let condition 14 no/a\ Tnen 



[— 7T,7r 



ld a(zd)F(zd,N e ,m d ) 
iz d 



x 



y- / c i Z2 i M z ^) F ( z ^ N ^ m ^ ( f ^h ^i^FjzuN^m!, 

|m 2 -f2|<2_ 7r |mi-^i|<2_ 7r 



xe ± ^ CT(e - 2)r/e n CT (0 - z)fl(e[r/e] - £miV e , 0)rLy (0)dziJdz 2 J . . . dz d )d9 + o r (l),(4.26) 
where o T (l) ^ as e — > /or any r G M 1 . 

Proof. We generalize the strategy of the proof of Proposition 3.6 from where this 
assertion is proved for d — 1 . The asymptotics (|4.26|) follows from ()4.25|) if we prove that 
the series over max^ \rrtj — Vj\ > 3 vanishes as e — > . 

First, let us consider the inner integral over z\ in (j4.25j) and denote it by I £ (6,z',m) : 

J e if+{e,z, mi )N e _ e if-{6,z,mx)N e 

I £ (9,z',m) = J a(9,z,m) : dz±, 

— 7T 

where a(9,z,m) = a(z 1 )g(9 — z)U a (9 — z)R(e[r/e] — emN £ , 9)H a >(9) , z' = (z 2 ,...,z d ) G 
[-7r,7r] d -\ 9 G [-n,n] d and 

f+(9,z,mi) = z 1 l l /N £ + z 1 (m 1 + 1/2) ±Lu a (9-z)T/(eN £ ), 
f-(9,z,7m) = z 1 l 1 /N £ + z 1 (m 1 -l/2)±co CT (9-z)r/(eN £ ). 

We have f + (9, 2, mi) =f_(9,z,m 1 ) = ±u a (9 1 , 9 2 - z 2 , . . . )t/(sN £ ) , and 



21=0 



21=0 



V 1 f ± {9,z,m 1 )\ zi=0 = l 1 /N e + m 1 ±l/2-u 1 . 
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Hence, Vi/±(#, z, mi) 7^ for \mx — Ui\ > 3. Indeed, we can admit that \h/N £ \ < 1 

21=0 

since N £ — > 00 as e — > , and the number Ii G Z is fixed. Further, we apply to / £ (0, 2', m) 
the limit of Lemma 3.7 from [2|, 



lim 



I £ (9, z', m) - 7re ± ^ (9l ' 92 - Z2 '- )r/£ a(0, (0, z'), m) 
x (sgnVi/+(0,z,mi)L =o ~ sgnVi/_(6», z, mi)| 21=0 



0. 



(4.27) 



Moreover, we obtain that I £ (8,z',m) — > as e — > , uniformly in 9 G [— 7r,7r] d and z 1 G 



since sgnVi/+(0, z, 



21=0 



21=0 



. We proceed by induction 



sgnVi/_(0,z,mi) 

for each inner integral over z 2 , ■ ■ ■ , Zd and obtain that the integrals with max^ \rrij — Vj\ > 3 
vanish as e — > . Further, we have to prove that the series over max.,- \rrij — Vj\ > 3 
also vanish. This follows from two facts: i) the function a(8,z,m) has a structure of 
f(9, z)R(e[r/e] — emN £ , 9) with a smooth function / , and ii) R(r, 9) satisfies condition 14. 
More exactly it is proved for the case d = 1 in [2] . The proof admits generalization to the 
case d > 1 , here we omit the detailed computations. 
Step v) The next step is to prove that 



(27T) 



-2d 



-i8-l+iui(9)T/e 



g{9)Ii a {9)R{r t Vuv(0)t, 9)U a ,(9) 



iz d l d +iv d N e z d 



atiZd) 



J5/2N E z d 



-i5/2N £ z d 



iz2h+iv2N E z 2 



J5/2N £ z 2 



-i5/2N £ z 2 



alz 2 



IZ 2 



:> iz\l\ J r iv\N e z\±iu) tT (9—z)T/£ 



aizi 



J5/2N e Z! _ ~i5/2N eZl 



IZ\ 



g{9 — z) dz\ ) dz% ) . . . dzd ) d6 



-oJl) 



(4.28) 



where o T (l) — > as e — > for any r G M 1 . It follows from ()4.26|) and the formula 
Y^\ m -v\<2 F ( z i' N ^ m i) = e iv * NeZ i{e ih l 2N * z i -e-^/^zj) . Formula KM is proved in Lemma 
3.8 from [2] for the case d = 1 . The proof is based on the condition 14 for function R and 
admits extension to the case d > 1 . 

Further, we apply ()4.27|) to the inner integrals from the RHS of (|4.28J) and obtain, for 
the inner integral over Z\ (denote it by I e (9, z') , where z' = (z 2 , . . . , z d ) ), 



lim \l F 



27ce ±i^(8i,e2-z2,..)T/s g ( du e 2 - z 2 ,... 1 e a 



0. 



(4.29) 



since in this case I £ (9,z') = f Tl a(zi)g(9 — z)(exp(if + (9,z)N £ )—exp(if_(9,z)N £ )^/(iz 1 )dzi 
with 

f + {9,z) = z 1 l 1 /N e + zi5/2 + u 1 z 1 ±w IT {0-z)T/(eN e ), 
/_(0,z) = z 1 l 1 /N £ -z l S/2 + p 1 z 1 ±uj a {9-z)T/{eN e ), 
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and sgnVi/±(0, z) = sgn(li/N e ±5/2) = ±1 for fixed l± E Z and small enough e > . 

21=0 

Finally, we obtain 

f = (2tt}_ / e -*v( w ^ w ^ w ) T/e 0(0)^ ( 430) 



4 



where o T (l) — > as e — > for any r ^0. 

Step vi) Note that the identities u a {9) ± av(#) = const ± in the exponent (see ()4.30jl ) 
with the const± 7^ are impossible by the condition E5. Furthermore, the oscillatory 
integrals with uj a {9) ± u a i(9) ^ const± vanish as e — > by the condition II and the 
Lebesgue-Riemann theorem. Hence, only the integrals with uj a {9) — u a i{9) = contribute 
to the integral ()4.30|) since ui a (9) + u a >(9) = would imply LU a (9) = to a '{9) = which is 
impossible by E4. We return to formula (|4.23|) and applying (|4.24|) one obtains formulas 
(P . ■ 



5 Proof of Theorems 13.41 and 13.5 



5.1 Convergence of Wigner matrices 

Proof of Theorem 13. 4L Theorem 14.11 implies that for any r E W 1 , r 7^ and y E (2Z) d 
the following convergence holds, 

lkE? /s (o([r/e]+i,/2)*®o([r/e]-i//2)) = W p (r;r,n), (5.1) 

e— >0 ' 

where in the Fourier space one has 

W p (r;r,fl) = I (^00^1/2 + n -i/2£i ^-1/2 

+in 1 / 2 q%.(9)n- 1 / 2 - in-^q^n 1 ^ 

= W p (r;r,9), (5.2) 

by formulas (l3~T2J) . (l3~T3l) and (pTTjl - ljOj) . Then convergence fEHij) follows from (pTT|l . (lo~2D 
and Lemma f5. 11 

Lemma 5.1 Lei conditions V2 and El - E3, E6 ao/d and a < —d/2 . Then 

sup sup \\Q £i t{x,y)\\ < C < 00. (5.3) 

Proof Applying (|4.18|) one has 

Q%(x, y) = W (Y\x, t) ® F% t)) = (Q £ (x\ y'), %{x\ t) ® t)), 

where 

^0',*) = (^(ar-afygffr-zO), ^ G Z d , i = 0, 1. 
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Then the Parseval identity, ()4.9|) and condition E6 imply 

\m-,t)\\l = (^r d j mo,t)\ 2 de = {2^y d J (\gf{9)\ 2 + \g\\e)\ 2 ) de<c < oo. 

Then Corollary 14.41 gives 

\Q%{x,y)\ = \(Qs(x',y'),K(x',t)®& y (y',t))\ < C\\*i(-M*W(',t)h<Ci < oo, 
where the constant C\ does not depend on x, y G Z d , t G R and £ > . ■ 

5.2 Weak convergence of measures [L e T i er as e — > 

Theorem 13.51 follows from Propositions 15.21 and 15.31 Proposition 15.21 ensures the existence 
of the limit measures of the family {{i £ T , r , e > 0} , while Proposition 15.31 provides the 
uniqueness. 

Proposition 5.2 Let conditions V2 and El - E3, E6 hold. Then for any r G M. d , r^O, 
the family of measures {^, £r , e > 0} is weakly compact in 7i a with any a < —d/2, and 
the following bound holds, 

sup / \\Y \\Iii £ t/ (dY ) < oo. (5.4) 

£>0 J 

Proof Definition (j2.1j) implies 

J \\Y4l^ /£ ^)=n{\\T^/e ] U{r/e)Y4l) 

= X> + \x\ 2 r(iiQZ /£ {[r/e] +x, [r/e]+x)+trQ% r/e {[r/e]+x, [r/e] +*)). 
xez d 

Since a < —d/2 , ()5.4|) follows from the bound (|5.3|) . Now the compactness of the measures 
family {fit, t £ K.} follows from the bound ()5.4|) by the Prokhorov Theorem Lemma 
II. 3.1] using the method of [18|. Theorem XII. 5. 2], since the embedding 7i a C Tip is compact 
if a > [3 . M 

Denote by Q T>r the quadratic form with the matrix kernel {(ff r {x — y))ij=o,i , 

Qr,r(*, = (^(ar-y),**^) ®* J '(y)), * G P. (5.5) 

«J=0,1 as,! 



Proposition 5.3 Let conditions VI - V4 and El - E6 hold. Then for any ret* 1 , 
and *6P, 

limy" exp(z(F, ^))^ /e>r {dY) = exp { - M>)} • (5.6) 

Proposition 15.31 is proved in Sections 6-9. 
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6 Convergence of characteristic functionals 

To prove Theorem 13. 5| it remains to check Proposition 15.31 Let us rewrite ()5.6)1 as 

#,.,,.(¥) = El(exp{ l (T^ r/£] U(r/e)Y ^)}) - e _> 0. (6.1) 

We will prove it in Sections 8, 9. In this section we evaluate (Ti r / e ]U(t)Y , ty) , t e R, by- 
using the following duality arguments. 

6.1 Duality arguments 

Remember that Yq G 7i a with a < —d/2 . For t G K introduce a 'formal adjoint' operator 
from space V to 7Y_ Q : 

(F, U'(t)V) = (U(t)Y, y G W a . (6.2) 

Let us denote by $ r (-,t) = U'(t)T[ r / E ]^f . Then using ()6.2|) we obtain 

(T- [r/e] U(t)Y ,*) = (Y ,$ r {;t)), i £l, £>0, rel d (6.3) 

The adjoint group U'(t) admits the following convenient description. Lemma 16.11 below 
displays that the action of group U'(t) coincides with the action of U(t) , up to the order 
of the components. 

Lemma 6.1 For \& = (\l/ , ^l/ 1 ) G V we have 

$(;t) = U'{t)* = (i>(;t)M-M ( 6 - 4 ) 
where if)(x,t) is the solution of Eqn \2.1\) with the initial data (uq,Vq) = (\E' 1 ,\1/ ) . 

The lemma allows us to construct the oscillatory integral representation for $> r (x, t) . 
Namely, ()6.4|) implies that in Fourier representation for $ r (-,i) = U'(t)T\ r / £ ]ty we have 

$ r (0,t) = A*(6)$ r (6,t), $(M) = g*(9)e t[r/£] - 6 ^(6), 

where 

( -V{6)\ cosS W -n(6) sin n(8)t\ 

^ {) \1 J' yt{) (6) sin Q{6)t cos tt(6)t )' 

Therefore 

$ r (x,t) = (2rr)- d J e- i6 - x g;{9)e i[r/e] - ^{9) d9, x G Z d . (6.5) 

T d 

Definition 6.2 D° = {$ e D : f («) = in a neighborhood of C} . 
From ()6.5|) we obtain 

Lemma 6.3 For any fixed \l/ G T>° the following bounds hold: 

i) \& r (x,t)\ < C r d / 2 , xeZ d . 

ii) For any p > there exist C p) ^ g > such that \$ r (x,t)\ < C p (\t\ + |x| + l)~ p , 
This lemma follows from Lemma 14.51 and the definition of T>° . 
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6.2 Equicontinuity of characteristic functionals 

Let us show that we can restrict ourselves to f £ P°. 

Lemma 6.4 The convergence h6. 1\) it suffices to prove for $ G D° only. 
Proof. Step i). For simplicity, let us put t = rje. Denote by 

2 J.,e 



sup |a A r(*,*)l<cK, m e v. 

£>o,teR,re 



We prove at first that 



Indeed, by (jfj.Hjl we have 

Qe,t,r(*, *) = Eg(|<r_ [r/e] E7(t)y , *> I') = (Qe(x, V ) , f) ® $ r (j/, t)>. 

So, by Corollary 14.41 we obtain 

sup |Q M , r (*,*)|<C sup ||$ r (-,t)||i. 

£>0,ieM,rS 



Finally, by the Parseval identity and condition E6, we get 



||*r(-,f)ll2 = (27r)- d / \\Gl 



i>{6)\ 2 d6 < C\\^\\l. 



Yd 



The bound ()6.6|) is proved. 

Step ii). By the Cauchy-Schwarz inequality, 



/$ P (*i)-/$ r (* 2 )| = \j "(e i(W -e^*))^ r (dr)|< J 



(6.6) 



e i(y,*i-* 2 > _ ]_ 



1/2 



= (Q £ ,t,r(*l - #2, ^1 - ^2)) V2 < C||*i - # 2 || 2 , 

where a constant C does not depend on e > , t 6 R 1 and r G M d . Hence, the charac- 
teristic functionals At/ er 0^)i ^ G M. , £ > , r 6 R 1 *, are equicontinuous in the space V 
endowed with the norm £ 2 . In the turn, the set V° is dense in this space. ■ 
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7 Bernstein's 'rooms-corridors' partition 

Let us introduce a 'room-corridor' partition of the ball {x G Z d : \x\ < j g t} with 7 g from 
Lemma Id 31 ii). For t > we choose below At, Pt G N (we will specify the asymptotical 
relations between t , A t and p 4 ) . Let us set h t = A t + p t and 

a* = jh, V = a? + A t , j G Z, n t = [7^/^]. (7.1) 

We call the slabs R\ = {x G Z d , |x| < : a J < < fr 7 } the 'rooms', C{ = {x G 

Z d , |x| < rif/ij : W < Xd < a J+1 } the 'corridors' and L t = {x G Z d , |x| > n t h t } the 'tails'. 
Here 2 = (x\, . . . ,Xd) , A t is the width of a room, and p t is that of a corridor. Let us 
denote by xl "the indicator of the room R\ , ^ that of the corridor C\ , and r\ t that of the 
tail L t . Then 

nt — l 

where the sum ^ . stands for . Hence we get the following Bernstein's type represen- 

j=-n t 

tat ion: 

(y ,<M-,*)> = y;.[(y ,^r(-,t)) + (n, £&■(•,*)>] + (r ,^$ r (-,t))). (7.2) 

Let us introduce the random variables r\ , , /f by 

ri = {Y ,xi$ r (;t)), 4 = (Y ,eM;t)), h = (Y , Vt $ r (; t)). (7.3) 
Then ()7.2|) becomes 

(Y ,$ r (.,t)) = J2.(rl + ci) + l t . (7.4) 

Lemma 7.1 Lei conditions VI - V2 ZioZd and \I/ G T>° . T/ie following bounds hold for 
t>l: 

K{\rl\ 2 ) < C(^) A t /t, Vj, (7.5) 

K(\c{\ 2 ) < C{H?)pt/t, Vj, (7.6) 

K{\h\ 2 ) < C P (V) (l + t)- p , Vp>0. (7.7) 

Proof The bound ()7.7j) follows from (j4.17|) . We discuss (|7.5|l . and ()7.6|) can be done in a 
similar way. Let us express Eg | 2 ) in the correlation matrices. Definition ()7.3)1 implies 
by the Fubini Theorem that 

E £ (|r?| 2 ) = (Q £ {x,y),xi{x)$ r (x,t) ®xi{y)®r{y,t)). (7.8) 
According to (ETTTj) and (ETTKl) . Eqn (TT51) implies that 

K(Wt\ 2 ) < cr d J2xKx)\\Q £ (x,y)\\ 

= Cr d J2xl(x)J2\\Qe(x,y)\\<CA t /t, (7.9) 

X z 

where ||(5 e (:c, y)\ stands for the norm of a matrix (Q % J(x, y)) . Therefore. ()7.9|) follows from 
Lemma 14.31 ■ 
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8 Ibragimov-Linnik Central Limit Theorem 

In this section we prove the convergence (j6.1|) . As was said, we use a version of the Central 
Limit Theorem developed by Ibragimov and Linnik. If Q Tr {^ 1 \l>) = 0, the convergence 
(1(1 lj) is obvious. Indeed, then, 

|E £ (exp{ ? (F ,$r(^r/£))})-^ r (vI/)| =E e (|exp{z(r ,^(-,r/ £ ))}-l|) 
< W (\(Y ^ r (;r/s))\) < (W (\(Y ,$ r (.,T/E))f)) 1/2 



((Q £ (x,y),^ r (x,T/e)®^ r (y,r/e))) l/2 = (Q £ , r/e , r (*, *)) 



1/2 



where Q e , T / e ,r(*, *) -> Qr,r(*, *) = , e -»• . Therefore, ffBTD follows from Theorem EHl 
Thus, we may assume that for a given \1/ e 22° , 



Q T , r (*,*)^0. 



Let us choose < 5 < 1 and 



(8.2) 



log* 



£ — ► oo. 



(8.3) 



Lemma 8.1 The following limit holds, 

'f) 1/2 + (i +ft H 



where a constant k > 



0, £ — > oo, 



U) 



Indeed, (18.31) implies that h t = p t + A t ~ , £ — ► oo . Therefore, n t ~ — ~ log £ . Then 

log£ h t 

dHUD follows by (Q. ■ 
For simplicity, we put t = t/e . By the triangle inequality, 



E^(exp{i(y ,^(-,£))})-A^W 



< 



E^(expO(y ,^(-,t))})-E^(exp{iJ].^}) 
ex P{ " ^£ /5(|ri| 2 )}-exp{ - ^Q r>r (*,*)} 
E £ ( exp{z^ ,r| }) -exp { - ^ Eg (K| 2 ) } 



= h + I 2 + h. 

We are going to show that all summands Ji , J2 , I3 tend to zero as £ — > 00 . 
Step (i) Eqn (|7.4|) implies 



5.5) 



J a = k(exp{z^.ri'}(exp{z^.c?' + ^}-l)) 



< 



1/2 



1/2 



(8.6) 
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From (JEHU), (ESJ), (J77|) and O we obtain that 

/i < Cn t (p t /t) 1/2 + C p r p -> 0, t -> oo. 
5Yep (mJ By the triangle inequality, 



i|^E £ (|r^| 2 ) - Q^ r (^,^)| < i |C A r(*,*) - Gr,r(*,*) 



J21 + I22 + -^23 5 



5.7) 



where Q e ,t,r is the quadratic form with the matrix kernel Q l £ J tr (x,y) . Theorem 13.41 implies 



that I 
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t=r/s 



as e — > . As for J 2 2 , we first obtain that 



/ 22 <E|E £ (r^)|- (8-8) 

The distance between the different rooms R\ is greater or equal to p t according to (J7TTJ). 
Then, by Lemma f6. 31 i) and condition V2, 

I22 < ^\(Q E {x,y)^ r {x,t)®^ r {y,t))\ 



3<l 



\-1 



j<i e y 
+00 

~ r d n 2 t t d - l A t J (1 + s^V" 1 rfs ~ + p t )" 7+d , 



(8.9) 



which vanishes as t — > 00 because of (|8.4|) and j > d. Finally, it remains to check that 
^23 — > , t — > 00 . We have 

Q £ , t>r (*, *) = Eg ((y , <&r (-, t)) 2 ) = K ( (E, (W + 4) + h) 2 ) , 

according to (|7.4|) . Therefore, by the Cauchy-Schwarz inequality, 

< |E e o((E^) 2 )- E o((E/*'+E^+^) 2 )| 

< Cr^E /5(l4 2 ) +^i(E £ ((E i ri') 2 )) 1/2 ( ri *E E o(l4i 2 ) +E £ (|Z t | 2 )) 1/2 
+CEg(|/ t | 2 ). (8.10) 

Then (JUJ), dSIE) and (JEHJ) imply 

E o((E/') 2 ) ^ E/od^'n+sE^J^H^) 

< CnAA + dn t (l + p0 _7+d < C 2 < 00. 
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Now dZD, (EH), fpnU|) and (JH3D yield 

/ 23 < Ci^pt/t + CW^/t) 1 / 2 + C 3 r p -> 0, t -> oo. 
So, the terms J21 , ^22 , ^23 in ()8.7|) tend to zero. Then ()8.7|) implies that for t — r/e 

^\ |E/o(Kf)-Q.,r(*,*) 

S^ep fmj It remains to verify that for t = r/e 
h 



0, £^0. 



(8.11) 



0, £ -> 0. 



Eg(«p{<X;.W'}) -exp{ - ^ E o(IW'| 2 )} 
Condition V3 yields 

n t — 1 

Eg(exp{i^7i'}) - J] E £ (expK}) 

nt— 1 

Eg(exp{ir t - ,H }exp{i ^ r£}) - E £ (exp{^})E £ (exp{ ? ^ rft) 
-jit+i 

nt— 1 n t — 1 

+ |E^(exp{zr-" t })Eg(exp{^ ^ rj}) - f] Eg(«p{irj}) 
< C(l + Pt )-* + 



< 



nt— 1 



-nt+1 



-nt+1 



nt— 1 



nt — 1 



E^(exp{z ^ r>'})- J] Eg(exp{irft) 

—nt+1 -nt+1 

We then apply condition V3 recursively and obtain, according to Lemma |8~T1 

nt— 1 



t=r/e 



0, £ -> 0. 



It remains to check that for t = t/s 

nt— 1 



J] E £ (exp{^}) -exp{ - ±£ E o(K| 2 )} 



0, £ -)• 0. 



According to the standard statement of the Central Limit Theorem (see, e.g., Theorem 
4.7]), it suffices to verify the Lindeberg condition: W5 > , 



CTt"— 'J 



0, £ -> 0. 



t=r/e 



Here a t = Xlj^o(l r t | 2 ) > an d ^o' a (/) = ^o(X a f) > where X a is the indicator of the event 
\f\ > a 2 . Note that (IRTTD and (Q imply that cx r/e -> Q^*,*) ^ 0, £ -> 0. Hence it 
remains to verify that 



Y^K'^K/f) - 0, £ - 0, for any a > 0. 
We check Eqn ()8.12|) in Section 9. This will complete the proof of Proposition 15.31 



(8.12) 
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9 The Lindeberg condition 

The proof of (|8.12J) is reduced to the proof of the following convergence 

E/o(K/J 4 )-0, 5-0, (9.1) 

by using Chebyshev's inequality. We deduce (j9.1j) from the following lemma. 

Lemma 9.1 Let the conditions of Theorem \3.£\ hold. Then for any ^ G T>° the following 
bounds hold, 

K(\rin<C(^)A 2 t /t 2 , t>l. (9.2) 

Proof. Step 1 Given four points G Z , we set 

Mf ) (x 1 , ...,x 4 ) = E £ (Yix 1 ) <g> ... ® Y(x A )) . Then, similarly to (JUJ) we have 

Wfl = (xKx 1 ) ■ ■ ■ xK^MPix 1 , . . . , z 4 ), $ r (x\ t) ® • • • ® $ r (x 4 , t)>. (9.3) 

Lemma f6.3| i) implies 

4 

E5(|rJ| 4 )<Cr M X; E ^(^•••^'(^ 4 )|Mi 4 )(x)i. (9.4) 

i=2 xe(z d ) 4 

By condition V4, we have 

xl^)... X l{x')\Mf\x)\< WW. (9-5) 

xe(z d ) 4 (n,*2,i3,*4)eP{i, 2,3,4} 

where 



Similarly to (|7.8jl . we have 

x E^^E**^ 1 + 1^ 3 -^ 4 i)~ 7 

~ [A^ 1 J^(l + \x h - x i2 |)" 7 2 . 

The sum in x 12 is bounded since 7 > d . Hence, 

VW.uW ^ C'At 2 ^" 2 - (9-6) 

Now the estimate (|9.2|) follows from (|9.4|) . ()9.5|) and (|9.6|) . This completes the proof of 
Lemma 19.11 ■ 
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